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Abstract 

For SU(3) lattice gauge theory we study properties of static quark 
sources represented by local Polyakov loops. We find that for tempera- 
tures both below and above T c coherent domains exist where the phases 
of the local loops have similar values in the vicinity of the center values 
0,±27r/3. The cluster properties of these domains are studied numeri- 
cally. We demonstrate that the deconfinement transition of SU(3) may 
be characterized by the percolation of suitably defined clusters. 
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Introductory remarks 

Confinement and the transition to a deconfining phase at high temperatures are 
important, but not yet sufficiently well understood properties of QCD. With the 
running and upcoming experiments at the RHIC, LHC and GSI facilities, it is 
important to also obtain a deeper theoretical understanding of the mechanisms 
that drive the various transitions in the QCD phase diagram. 

An influential idea is the Svetitsky-Jaffe conjecture [1J which states that 
for pure gluodynamics the critical behavior can be described by an effective 
spin model in 3 dimensions which is invariant under the center group (for 
gauge group SU(3)). The spin degrees of freedom are related [2J to static quark 
sources represented by Polyakov loops, which in a lattice regularization are given 
by 

N 

L{x) = tr c JJ Ckfo *) • (!) 
t=i 

The Polyakov loop L{x) is defined as the ordered product of the SU(3) valued 
temporal gauge variables U±(x,t) at a fixed spatial position x, where N is the 
number of lattice points in time direction and tr c is the trace over color indices. 
The loop L(x) thus is a gauge transporter that closes around compactified time. 
Often also the spatially averaged loop P = l/V^gL(x) is considered, where 
V is the spatial volume. Due to translational invariance P and L(x) have the 
same vacuum expectation value. 

The Polyakov loop corresponds to a static quark source and its vacuum 
expectation value is (after a suitable renormalization) related to the free energy 
F q of a single quark, (P) oc exp(— F q /T), where T is the temperature (the 
Boltzmann constant is set to 1 in our units). Below the critical temperature T c 
quarks are confined and F q is infinite, implying (P) = 0. This is evident in the 
Ihs. plot of Fig.[TJwhere we show scatter plots of the values of the Polyakov loop 
P in the complex plane for 100 configurations below (Ihs. panel) and above T c 
(rhs.) 0- In the high temperature phase quarks become deconfined leading to a 
finite F q which gives rise to a non-vanishing Polyakov loop (rhs. in Fig. [TJ. 

On a finite lattice, above T c the phase of the Polyakov loop assumes values 
near the center phases which for SU(3) are 0, ±2-7r/3 (rhs. plot of Fig. [T|. This is 
a reflection of the underlying center symmetry which is a symmetry of the action 

1 The numerical results we show are from a Monte Carlo simulation of SU(3) lattice 
gauge theory using the Liischer-Weisz gauge action [3]. We work on various lattice sizes 
ranging from 20 3 x 6 to 40 3 x 12. The scale was set [4] using the Sommer parameter. 
In our figures we always use the dimensionless ratio T/T c with the critical temperature 
T c = 296 MeV calculated for this action in 5]. All errors we show are statistical errors 
determined with a single elimination jackknife analysis. 



1 





1 1 


1 


1 i 1 i 

T = 0.63T - 

c 


1 1 
- 


1 1 1 1 1 
T=1.32T - 

c 


0.2 












0.0 




o 








0.2 








- 9 




0.4 


I . 


1 


, 1 , 1 


i . 


I , I , I 


-0.2 


0.0 


0.2 0.4 


-0.2 


0.0 0.2 0.4 



0.2 
0.0 
-0.2 
-0.4 



Figure 1: Scatter plots of the spatially averaged Polyakov loop P in the 
complex plane for configurations below (lhs. panel) and above T c (rhs.). We 
show the results for our 40 3 x 6 ensembles. 



and the path integral measure of gluodynamics, that is broken spontaneously 
above the deconfinement temperature T c . As long as the volume is finite all 
three sectors are populated, while in an infinite volume only one of the three 
phase values survives. This center symmetry and its spontaneous breaking are 
the basis for the above mentioned Svetitsky-Jaffe conjecture [TJ. 

The relation of the deconfinement transition of SU(iV) gauge theory to Zjy- 
symmetric spin models has an interesting implication: For such spin models it 
is known that suitably defined clusters made from neighboring spins that point 
in the same direction show the onset of percolation at the same temperature 
where the Zjv-symmetry is broken spontaneously. For, e.g., the Ising system 
these percolating clusters were identified [6] as the Fortuin-Kasteleyn clusters 
|7j. An interesting question is whether the cluster- and percolation properties 
can be directly observed in a lattice simulation of gluodynamics - without the 
intermediate step of the effective spin theory |2J for the Polyakov loops. 

For the case of gauge group SU(2) such cluster structures were analyzed in 
a series of papers [8j [9], while for SU(3) the formation of center clusters has 
not yet been explored. In this paper we try to close this gap and study the 
behavior of the local loops L{x) and the possible formation of center clusters. 
Furthermore, we study center clusters not only near T c (where they directly can 
be expected from the Svetitsky-Yaffe conjecture) but explore their emergence 
and properties in a window of temperatures ranging from 0.63 T c to 1.32 T c . 
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Figure 2: Histograms for the distribution of the modulus p(x) of the local 
loops L(x) for temperatures below and above T c . The full curve is the 
distribution according to Haar measure (40 x 6 ensembles). 

Properties of local Polyakov loops 

For analyzing spatial structures of L{x) on individual configurations we write 
the local loops in terms of a modulus p(x) and a phase <p(x), 

L(x) = p(x)e 1 ^ . (2) 

The first step of our investigation is to study the behavior of the modulus p(x). 
In Fig. [2] we show histograms for the distribution of p(x) in the confined (Ihs. 
plot) and the deconfined phase (rhs.). It is obvious, that the distributions of the 
modulus p(x) below and above T c are almost indistinguishable. Furthermore we 
find that the distribution follows very closely the distribution according to Haar 
measure, which we show as a full curve. Only above T c we observe a very small 
deviation from the Haar measure distribution. The Haar measure distribution 
curves for the modulus and the phase are defined as 

P(p) = J D[U]6(p-\Tr[U]\) , P(ip) = J D[U] 8(<p — arg Tr [U]) , (3) 

where 6 is the Dirac delta-function and D[U] is the Haar integration measure 
for group elements U G SU(3). These two distributions are obtained from a 
single group element and thus do not depend on any lattice parameters. 

From the fact that the change of the modulus is very small we conclude 
that the jump of (P) at T c , signaling the first order deconfinement transition, is 
not driven by a changing modulus of the local loops L{x). Thus we focus on 
the behavior of the phase (p(x), and again study histograms for its distribution. 



3 



H[q>(x)] 
0.6 

0.4 

0.2 





l 1 1 1 1 l 1 1 1 1 l 1 1 1 1 l 1 1 1 1 l 1 1 1 1 l 1 1 1 1 l 



T = 0.63 T 




i .... i .... i .... i .... i .... i .... i 



i 1 1 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i 1 1 1 1 i 1 1 1 1 1 



T = 1.32 T 



lMltlTTTTTTJ» 



i .... i .... i .... i .... i .... i .... i 



H[<p(x)] 

0.6 

0.4 
0.2 




-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 

cp(x) 



-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 

q>(x) 



Figure 3: Histograms for the distribution of the phase tp(x) of the local loops 
L(x). We compare the distribution below T c (lhs. plot) to the distribution 
in the deconfined phase (rhs.) for the sector of configurations with phases 
of the averaged loop P near —2ir/3. The full curve in the lhs. plot is the 
distribution according to Haar measure (40 x 6 ensembles). 



In Fig. [3] we compare the distribution below T c (lhs. plot) to the one in the 
deconfined phase (rhs.). For the latter we show the distribution for the sector 
of configurations characterized by phases of the averaged Polyakov loop P in 
the vicinity of —2ir/3 (compare the rhs. of Fig. [1]). 

The distribution of the phases <p(x) is rather interesting: Also in the confined 
low temperature phase (lhs. plot in Fig. [3} the distribution clearly is peaked at 
the center phases — 2tt/3, and +2-7T/3, and again perfectly follows the Haar 
measure distribution (full curve in the lhs. plot). The distribution is identical 
around these three phases and the vanishing result for (P) below T c comes from 
a phase average, 1 + e i27r / 3 + e _i27r / 3 = 0. 

Above T c (rhs. plot in Fig. [3]) the distribution singles out one of the phases. 
In our case, where configurations in the sector with phases of the averaged loop 
P near —2ir/3 are used for the plot, it is the value — 2tt/3 which is singled out. 
For configurations in one of the other two sectors (see rhs. plot of Fig. [I]) the 
distribution is shifted periodically by ±2-7r/3. Obviously, above T c the distribu- 
tion is not equal for the three center phases and the cancellation of phases does 
no longer work, resulting in a non-zero (P). 

The histograms for the phases tp(x) suggest that at the critical temperature 
the local loops L(x) start to favor phases near one spontaneously selected center 
value, while phases near the other two center values are depleted. This is 
illustrated in more detail in Fig. [4] where we show the abundance A of lattice 
points with phases of L(x) near the dominant and subdominant center values. 
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Figure 4: Abundance of lattice points x with phases of the local loop L(x) in 
the dominant (triangles) and subdominant center sectors (circles, squares) 
as a function of temperature (40 x 6 ensembles). 



To define the abundance A we cut the interval (— ir, n) at the minima of the 
distribution of Fig. [3] into the three sub-intervals (— it,— 7r/3), (— 7r/3, tt/3), 
(7r/3,7r), which we refer to as "center sectors" . We count the number of lattice 
points with phases in each of the three center sectors and obtain their abundance 
A by normalizing these counts with the volume. Fig. 0] shows that at low 
temperatures all three center sectors are populated with probability 1/3. Near 
T c one of the sectors starts to dominate while the other sectors are depleted. 

Coherent center domains 

We have demonstrated for a wide range of temperatures that the center sectors 
play an important role for the phases ip(x) of the local loops L(x), which 
cluster near the center phases 0, ±27r/3 at all temperatures. The deconfinement 
transition is manifest in the onset of a dominance of one spontaneously selected 
center sector. We now analyze whether the values of the phases tp(x) are 
distributed homogeneously in space, or if instead there exist spatial domains 
with coherent phase values in the same sector. 

In order to study such domains, we use sub-intervals that divide the interval 
{—it, it) for the values of the ip(x). For a more general analysis we introduce the 
cutting parameter 5 > and define the three sub-intervals as (— 7T+5, —it/3— 5), 
(— 7r/3 + S, 7r/3 — 5) and (tt/3 + 5, it — S) (which we again refer to as "center 
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Figure 5: Weight of the largest center cluster as a function of temperature. 

sectors"). For 5 = we obtain the old sub-intervals, while a value of 5 > 
allows us to cut out lattice points where the phases are near the minima of the 
distributions shown in Fig. [3] 

The definition of the clusters slightly differs from those that have been used 
for the analysis [8l[9] of clusters in SU(2) gauge theory. Besides modifications of 
the Fortuin-Kasteleyn prescription studied in [8J, in |9J the bonding probability 
between neighboring sites with same sign Polyakov ioop£| was introduced as 
a free parameter. This parameter could then be tuned such that the onset of 
percolation agrees with the deconfinement temperature. In our definition the 
parameter 5 allows one to reduce the lattice to a skeleton of points with phases 
close to the center elements (in intervals of width 2(tt/3 — 5) around the center 
values). For the plots shown in Figs. [5] and [6] we choose 5 such that roughly 
those 40 % of lattice points are cut where the phases do not strongly lean 
towards one of the center values. We found that near T c the critical properties 
of the clusters (behavior of largest cluster and percolation) are stable when S 
is varied in a small interval around that value [llj (compare also the discussion 
in |9J). For example the curve for the weight of the largest cluster (see Fig. [5] 
below, where we show a comparison of different spatial volumes for a cut of 39 
%) is form-invariant in a range of cuts from 30 % to 45 % and only is rescaled 
by a change of the amplitude of less than 15 %. 

In a next step we define clusters by assigning neighboring lattice sites with 

2 For SU(2) the L(x) are real and the center phases are either +1 or —1. 
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Figure 6: Percolation probability of the dominant center clusters as a 
function of temperature. 

phases ip(x) in the same center sector to the same cluster. Once these center 
clusters are defined we can study their properties and behavior with temperature 
using concepts developed for the percolation problem [10]. In Fig. [5] we show 
the weight (i.e., the number of sites) of the largest cluster as a function of 
the temperature. For low temperatures all clusters are small, while as T is 
increased towards the deconfinement temperature the largest cluster starts to 
grow quickly and above T c scales with the volume. This property indicates that 
in the deconfined phase the system has developed a percolating cluster. The 
onset of percolation at T c is confirmed in Fig. [6]where we show the percolation 
probability poo as a function of T/T c . The percolation probability is computed by 
averaging an observable which is 1 if a spanning cluster exists and otherwise. 
In our case, where we have periodic spatial boundary conditions, a spanning 
cluster is defined as a cluster who has at least one member site in every y-z 
plane. In other words, we analyze percolation in x direction, which is, however, 
no loss of generality as we have invariance under discrete spatial rotations. 
Varying 5 in a range where the number of points we cut varies between 30 and 
45 % (Fig. [6] is for 39 %), slightly roundens the transition curve, but leaves the 
onset of percolation unchanged at T/T c = 1. 

An interesting question is the size of the clusters in physical units in the 
confining phase, which could be related to some hadronic scale (see also the 
discussion in the next section). In order to study this cluster size below T c we 
computed 2-point correlation functions of points within the individual clusters. 
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These correlators decay exponentially oc exp(— r/£) with distance r, and the 
factor £ defines a linear size of the clusters in lattice units. We then analyze 
d = 2£a, which gives a definition of the cluster diameter in physical units (a 
is the lattice spacing in fm). We find that up to T = 0.85T C this diameter is 
essentially independent of the temperature, with a value of d = 0.46(5) fm at 
a cut of 39 % and d = 0.62(7) fm at a cut of 30 %. Compared to the expected 
sizes of rougly 0.5 fm for heavy quark mesons this is a quite reasonable result 
for the linear scale of the clusters which suggests that the physical role of the 
clusters below T c should be studied in more detail (see [11J). 

Summary and discussion 

We have explored the clustering of the phases ip(x) of the local quark sources 
L(x) near the center values, both below and above T c . We find that in the range 
of temperatures we consider, T = 0.63 T c to T = 1.32 T c , the local Polyakov 
loop phases prefer values near the center values and corresponding clusters may 
be identified for these temperatures. Using the parameter 5 we can construct 
clusters such that the deconfinement transition is characterized by percolation 
of the clusters in the dominant sector. 

From the cluster properties a simple qualitative picture for confinement and 
the deconfinement transition emerges. Below T c the clusters of lattice points 
which have the same center phase information are small. Only if a quark- and 
an anti-quark source are sufficiently close to each other they fit into the same 
cluster and can have a non-vanishing expectation value. Sources at distances 
larger than a typical cluster size receive the independent center fluctuations from 
different clusters and the correlator averages to zero. Above T c the clusters 
percolate and coherent center information is available also for larger distances 
allowing for non-vanishing correlation at large separation of the sources. In this 
picture deconfinement is a direct consequence of a percolating center cluster. 

A possible role of local center structures for confinement has been addressed 
also in a different approach, using a projection of the link variables U^(x,t) at 
all points in space and time to a center element after fixing to a suitable gauge 
(see, e.g., [12] for a selection of recent results). This analysis is motivated by 
understanding the role of topological objects for the QCD phase transition. It 
would be highly interesting to study a possible connection of the percolation 
aspects of the transition to the dynamics of such topological objects. Of par- 
ticular relavance would be an analysis of a possible relation to calorons which 
induce strong local variations of the Polyakov loop that might play an important 
role in the formation of the center clusters [13] . 

We conclude with a few comments on the extension of the center domain 
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picture to the case of full QCD: The fermion determinant describing the dy- 
namical quarks can be expressed as a sum over closed loops, which may be 
viewed as generalized Polyakov loops and are sensitive to the center properties 
of the gauge fields [14]. The fermion determinant breaks the center symmetry 
explicitly and acts like an external magnetic field which favors the real sector 
(phase 0) for the Polyakov loop P. However, preliminary numerical results with 
dynamical fermions pj] show that locally also the two complex sectors (phases 
±27r/3) remain populated. The corresponding clusters will again lead to a 
coherent phase information for sufficiently close quark lines. As for the pure 
gauge theory studied in this letter, the preliminary results [11] show that the 
transition to confinement is again accompanied by a pronounced increase of the 
abundance for the dominant (i.e., real) sector. However, the explicit symmetry 
breaking through the determinant leads to a crossover type of behavior in the 
dynamical case. An interesting related question, which has already been raised 
in the literature [15 j . is whether also the chiral transition may be characterized 
as a percolation phenomenon. 
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